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Abstract 



In this paper we give an algorithm for finding the mth base-6 digit 
of a positive integer n (m = 1 is the least significant digit) defined 
as the final number in a sequence of integers gotten by multiply- 
ing, adding, and subtracting previous numbers in the sequence (ac- 
tually, the algorithm finds arbitrarily precise approximations to n/b"^ 
(mod I), which can be used to get this mth digit whenever the lower 
m — 1 digits do not begin with a long run of the digit 6 — 1). In 
many cases, this algorithm will require far less memory than it takes 
to write down the base-6 digits of n, and will have a running time (in 
bit operations) only slightly worse than linear in the number of these 
base-6 digits. 

One easy-to-state consequence of the above result is that the mth 
base- 10 digit of 2 can be found using 

0(t2/3log^t) bits of memory 
and O(tlog^t) bit operations, where C > is constant. Of course, 
if m = 0(t2/3) then one can do much better (by just computing n 
(mod 6"*)), so the result is only non-trivial when t = 0{m?/'^). 

The algorithm we give is highly parallelizable, although if one uses 
M processors, to get an M-fold reduction in running time, the memory 
requirements will increase by a factor of M. 

1 Introduction 

Suppose that a is a positive real number. Then, a standard fact regarding 
base-6 representations is that given any integer b > 2, there exists and integer 
J, and a sequence of integers rj,rj_i, such that 




(1) 



1 



and we write 

a = (rjrj„i...ro.r_i...)fe 

to denote this expansion (when J > 0). If we further disallow r„ = 6 — 1 for 
all n < N, for some N, then the sequence of rj's is unique. If we make no 
such restrictions on r„, then the sequence of rj's may not be unique, as in 
the base- 10 expansion 

1 = 0.9999999... 

A question which has attracted recent attention (see (Jj) is the follow- 
ing: Given a base b, and an integer m > 0, efficiently determine a good 
approximation to 

u = {6™7r}, 

where for a real number 6, the notation {6} means 

{0} = e - ie\ (Note : {e} = e (mod 1), < {^} < 1). 

We will say that such an approximation 7 is a level-p approximation if and 
only if 

1 

' ' bP 

Now, if the approximation 7 is sufficiently good, then one can use it to 
determine the digits of vr: For example, if we take m = 3 and b = 10, then 
we get 

u = b'^n = lO^TT = 0.59265358979... (mod 1). 
Suppose we had an approximation 7 to satisfying 

1 

n — ul < : 

" ' 1000' 

so, 7 = 0.59.... Then, the leading base-10 digit of 7 is the same as the fifth 
base- 10 digit of tt (from the left), which is 5. 

More generally, a good approximation to u gives us the (m + 2)nd digit 
of vr; however, depending on the value of m selected, this approximation may 
need to be extremely close to u, in order to determine this digit. For example, 
in the above instance with m = 3, if 7 = 0.6, then 



|7 - '^l < .^o' 



1 

Too' 



2 



and we note that the leading digit of 7 is not the same as the leading digit 
of V. 



In this paper we will describe a method for determining the mth base-6 
digit of an integer, where m = 1 corresponds to the least significant digit 
(note that for the digits-of-vr problem above, m = 1 corresponded to the 
leading digit; so, the mth digit is defined differently in this context). As in 
the problem concerning digits of vr, this method produces arbitrarily precise 
approximations to 



where n is this integer, which will be defined by a certain type of expansion 
that we will describe below. Note that the value of m is at least 1 (if m < 0, 
then u is trivially 0). Now, if 



then the mth digit of n equals r_i; and, if A is a sufficiently good approxi- 
mation to z/, then the leading digit of A will also equal r_i. 

The type of expansion for n we will use is defined as follows: In j2], 
Smale and Shub say that a computation of length L of a positive integer n is 
a sequence of integers si, S2, sl, where si = 0, S2 = 1, and for i >3, 



where o is either addition, subtraction, or multiplication, and where sl = n. 
If one expresses such a computation as a string indices (j, k) and operations 
+, — and X, then given an integer N > 2, one can compute n (mod A^) 
using only 0(L(logL)(log A^)(loglog A^)^) bit operations, by just computing 
the sequence si, ...,sl modulo A^. The factor (log A^) (log log A^)^ in this big- 
O appears because a product of integers modulo A^ can be computed using 
Fast Fourier Transforms using only 0((log A^)(loglog A^)^) bit operations. 

We will say that an integer n has computational complexity L if and only 
if there exists a computation of lenth L for computing n. 

In the next section we will prove a general result, which we will use to 
prove the following theorem: 

Theorem 1 Given a positive integer n having computational complexity L, 
an approximation A to the number of base-b digits of n (see the input specs 




V = 



(0.r_ir_2 ■ ■ ■)b, 



Sj o Sfc, where j, k < i, 
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below), an integer m > 0, and a level y, there exists an algorithm for com- 
puting a level-y approximation to 



V 



n 1 



This algorithm requires only 

0{yL{log^^^ n) \og^{b + m + y + L + logn)) bits of memory, 

and 

0{yL{\ogn) \og^{b -\- m -\- y -\- L -\- logn)) bit operations. 

The input and output requirements of this algorithm are as follows: 

Input: The integers m, y, and a string representing the length-L 
computation ofn. Also, the algorithm requires as input an integer 
A, which is an approximation to the number d of base-b digits of 
n. This approximation need only satisfy 

1 ^ 

< -r < 2. 

2 d 

We further restrict m so that m < 2 A + y + 1, since otherwise 
7 = satisfies the conclusion of our Theorem. 

Output: The level-y approximation to v, encoded as a string of 
y + 1 base-b digits. 

This theorem requires a httlc more explanation. First of all, the input 
to the algorithm will be a string of 0(LlogL + logm + logy) bits, which 
is smaller than the space requirement listed above (when C > 1). The 
L log L, log m, and log y terms here account for the number of bits needed 
to specify the length L computation of n, the index m, and the level y, 
respectively. The approximation A to d requires only O (log logn) bits of 
space, and this turns out to be O(LlogL): To see why this is so, we note 
that any length L computation produces an integer n < 2^ , which can be 
proved by induction. It follows then that log logn = 0{L). 

The output of the algorithm will be a string of 0{y) bits, representing 
the base-6 approximation 7 to u. The number 7 will have only y base-6 
digits, and will satisfy 

1 

" ' by 
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Perhaps the most surprising aspect of the above theorem is that the indi- 
cated algorithm can require significantly less memory to find the approxima- 
tion to V than it does to write down the number n, which will have 0(log?7,) 
base-6 digits. The following corollary of the above theorem will make this 
point clear: 

Corollary 1 Suppose that a,b > 2 and m,t,y > 1 are all integers. There 
exists an algorithm which computes a level-y base-b approximation 7 to u, 
where 

a* 



This algorithm requires only 

0{yt^^^ log*" (a + b + y + m + t)) bits of memory, 

where C > 0, and performs 

0{yt log*" (a + b + y + m + t)) bit operations. 

Now, the number of bits needed to write down the number tt, = a* is clearly 
0{t log a); and yet, if, say, we take y = I, this algorithm requires only ^2/3+0(1) 
bits of space. 

This corollary follows since a* has computational complexity 

L = 0(\og\a + t)). 

To see this, we note that a* can be generated by repeated squaring: If t = 
2*^ + ■ ■ ■ + 2*" , then a* = ^ ■ ■ ■ " . These numbers can be computed by 
starting with a; then squaring to get a^; then squaring again to get a^; then 
continuing, this produces the list after only h multiplications. 



The rest of this paper is organized as follows: In the next section we will 
state the Main Theorem (Theorem |2)) and then use it to deduce Theorem 
n In section El we give a proof of the Main Theorem. Finally, in section 0] 
we give a proof of a proposition (Proposition^, which is an auxiliary result 
needed for the proof of the Main Theorem. 
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2 Main Theorem and Proof of Theorem [T] 



Theorem ^ is actually a corollary of a more general result concerning approx- 
imations to V. In this section we will state this result, which will henceforth 
be called the Main Theorem, and then show how to apply it to prove Theo- 
rem n The proof of this theorem, as well as a brief description of the ideas 
used to prove it, can be found in sectional also, at the end of Section ITT] 
we will give a brief statement on how to parallelize the algorithm. 

We suppose that 6 > 2 is an integer, which is to be the base used; that 
< ri/a* < 1 is some rational number where n, a > 1, t > are integers, 
and where n has computational complexity L; that ?/ > 1 is some level of 
precision to be used; and finally, that /i > is some integer. Then, given any 
pair of integers 5*, T satisfying 

ST > ^i^o^r^ + i^. + y + 2)logb)^ ^ 
logo 

we have the following 

Theorem 2 (Main Theorem) Let 

-=W^}- (3) 



There exists an algorithm which computes a level-y approximation 7 to this 
number v, where the space and time requirements of the algorithm are as 
follows: 

Space: 0{yL{S + T) \og^ {y + L + S + T + a + h + n + \ogt)) hits 
of memory, where C > 0. 

Time: 0{yL{ST + T^) \og^{y + L + S + T + a + b + fi + hgt)) 
bit operations. 

Here we give more precise information about the input and output spec- 
ifications of the algorithm: 

Input: The positive integers a, t, y, S and T, as well as a string 
of O(LlogL) characters representing the length-L computation 
needed to produce n. 
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Output: The algorithm will give an approximation 7 to u. This 
approximation will have y + 1 base-6 digits, and will satisfy 



I7 — z/| < 



by' 



To prove Theorem ^ using this result, we let a 
We note that this gives 

h^n n 

We also let 



6, it = 2^4, and /z = t—m. 



S 



T 



3(log&)(3A + /i + 7/ + 2) ^| 
log a 



2/3 



3(log&)(3A + /i + y + 2) 
log a 



1/3 



+ 1, 



+ 1. 



We note that this choice of S and T satisfies (j2I). ^ 

Now, applying the algorithm described in Theorem |21 with the parameters 
indicated above, we get the same output as described in Theorem [TJ The 
running time and space requirements to run this algorithm are also as stated 
in Theorem^ for our particular choices of S and T. 



3 Proof of Theorem |21 

Let S and T be as in (j21), and let r and k be integers such that 

t = Sk-r, 0<r<S-l. 

Then, we have that 

n noJ' 

The idea of the proof of Theorem |21 is to approximate h^a (and therefore 
z/) as follows: 

h^'a = 7i + --- + 7t + ^, (4) 

^To show this, one needs the fact that 3Alogb > logrt, which follows since n has < 2 A 
base-& digits. 
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where 

7,- = 



- J 



for some rationals Ai, At, and where 

\E\ < for ST sufficiently large. (5) 

Then, we will find approximations 7j,...,7^ to {71}, {7t}. Now, if the 
precision of these approximations 7^^, ...,7^ is high enough, and if we we let 
E satisfy 

S = {7; + 72 + --- + 7t}> 

then E will be an approximation to u; and, if we then take 7 to be the closest 
number to E having y + 1 base-6 digits, then 7 will be a level-?/ approximation 
to z/. 

We claim that the approximations 7j^,...,7^ to {71}, {7^} need only 
have 

loffT 



w = y + 



+ 3, 



log 6 

base-6 digits (and be level-w approximations), in order to guarantee that E 
is a level-y + 1 approximation to u. Note that this would imply that 

To see only w base-6 digits are needed, we note that if these numbers 7^ 
satisfy this last inequality, then by the triangle inequality. 



T 

I7-H < ^ + < ^ + + 1^1 

1111 

as claimed. 

Let us now find a set of values for Ai, At which make hold: Using 
the geometric series identity, we have that 



A- °° 5- 

1=1 j=i 
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where 



+ A22^-' + A^y-' + ■■■ + AtT^-\ 



We seek values for Ai, At so that 



0, 



for 1 < j < T, j T^k] 



na'^, for j = k. 
The following Proposition gives the solution we seek 
Proposition 1 We have that 



A, = 

and, 

Note that this implies 



na^ ^Coef. of ^ in ]^ 



h=i X — ri 



]\h=i J - h 



SaT 



\Aj\ < nTa^'i 



7i 



6^- 



-l)^-%a'' fCoef. of x'''^ in Ff^^i 



--1 X — h 



(j-l)!(r-j)!(a^-j) 
From this proposition we deduce that 

Ai + A22^'^ + ■■■ + AtT^~^ 



\E\ 



j=T+l 



< 



{nTaH^)T^ 



j=T+l 

aS{T+i) 1^ ^sj 

j=0 



(7) 



(8) 



for ST large enough; and so, (0) follows. 

We now have all the ingredients necessary to prove Theorem |2l which we 
will give as the following algorithm: 
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3.1 Algorithm 1 



The input, output, and requirements of this algorithm are as stated in The- 
orem |21 Here are the steps of the algorithm: 



1. Let 



w = y + 



logT 
log 6 



+ 3. 



Note that this choice of w satisfies 

T 



< 



2. Set S = 0, and let r, 5, T and k be as described at the beginning of 
this section. 

3. For j from 1 to T do steps 4 through 8. 

4. Compute 



5. Set 

V ^ {j-l)\{T-j)\Q. 

This number can be computed using 0{{T + S) log^(T + S)) bit operations, 
and just as much memory (for some C > 0). 

6. Apply Algorithm 2 (given in the next subsection of the paper) to 
compute 



H ^ Coef . of x^-^ in fj x - /i. 



This step requires 0(T^ log*" T) bit operations and 0(T log*" T) bits of mem- 
ory. 

7. Compute 

u <— i—l)^~^naJ' H (mod v), < m < i; — 1. 
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Since n has computational complexity L, this step requires only 0{L(S + 
T) log*" (L -\- S -\- T -\- a)) bit operations, and just as much memory. 



8. Find a number r having w base-6 digits satisfying 



where 



|r-0| < 



(So, T will be a level-ty approximation to (f).) 

We note that this number r can be easily computed by first letting 



Then, by finding the first w + 1 significant digits of uq/v^ and letting r = 
(0.r_i...r_.u,_i)b one see that the above inequalities are satisfied. 



We only need to do level-w arithmetic in base-6 here. 
(If j < T, then increment j and loop back to step 4.) 

10. (We assume j = T .) Let 7 be the number having y + 1 base-6 digits 
which comes nearest to E, and then OUTPUT 7. 

We note that we can perform the operations in steps 3 through 8, with 
different values of j, in parallel. For example, given two processors, we can 
assign processor 1 to perform steps 3 through 8, with values of j < T/2, and 
then assign processor 2 to do the same, but with T/2 < j < T. This would 
result in an two-fold reduction in the running time, as long as is sufficiently 
large. Of course, the memory requirements would double, because each of 
the two processors would require their own seperate memories. 

More generally, we have that, given M processors, for n sufficiently large. 
Algorithm 1 can be computed in parallel, resulting in an M-fold reduction 
in running time, but an M-fold increase in memory requirements. 



(mod v), < Uo < V — 1, 



and then noting that 




9. Set 



{S + r}. 
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3.2 Algorithm 2 



Input: T,k,j. 

Output: Coef. of x^~^ in n^=i x — h. 

Requirements: The algorithm performs O(T^log^T) bit opera- 
tions (for some C > 0), but requires only 0(Tlog*^T) bits of 
memory. 

1. Let P be the least integer such that 

p<P 
p prime 

Note: P — O(TlogT), and can be computed using O(Tlog^T) bit opera- 
tions (for some D > 0); and so, we can compute and store P within the time 
and space requirements listed above for the algorithm. We also note that 
every coefficient of the polynomial in the output specifications is less than 
A/2 in absolute value. 

2. Set S = 0. 

3. For each prime p < P do steps 4 through 8. 

4. Compute the polynomial 

T 

f{x) = JJ^ a; — /i (modp). 

h=l 

Note: This polynomial can be stored as a length-T coefficient vector, and 
the number of bits required to store such a vector is 0(T log p) = 0(T log T); 
also, this polynomial can be computed using O(Tlog^T) bit operations by 
making use of FFT's and a divide-and-conquer strategy for polynomial mul- 
tiplication. The divide-and-conquer part of the algorithm can probably best 
be described as the following recursive procedure: First, we suppose that L 
is a set of polynomials to be producted together modulo p, and Product (L) 
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denotes the procedure for computing this product. The pseudocode for this 
procedure is given as follows: 



If \L\ — 1 (i.e. L has only one polynomial), then 

RETURN the contents of L (mod p) ; 
Else, if |L| > 2, say L = {/i, ft}, then 

RETURN 

Product({/i, /Lt/2j})-Product({/Lt/2j+i, ft}) (mod p) 

Now, using FFT's to perform the polynomial multiphcation in this second 
step (the 'Else' step), we see that if the two polynomials being multiplied 
together have degrees £i and £2, respectively, then the multiplication should 
take no more than 0( (£1 + ^*2) log^(£i + £2+p)) bit operations. Now, if we run 
Product (L) starting with L consisting of all linear factors x — h, 1 < h <T, 
h j, then if T — 1 is a power of 2, the procedure products together (T — 1)/2 
pairs of degree 1 polynomials; (T — l)/4 pairs of degree 2 polynomials; and 
so on, all the way down to two polynomials of degree (T— l)/2. So, the total 
number of bit operations required to run this producedure is 



The memory requirements (in bits) are likewise of the same order. 



« (log^T) Yl 



-2^ = O(Tlog^r). 



T 



J<(logT)/log2+l 



5. Set 



H 



Coef. of x'^ ^ in f{x) (mod p). 



6. Set 



N 



(A/p) H (mod p), where < < p — 1. 



7. Set 



NA 



E 



E + 



P 
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8. Increment the value of p, and return to step 4, unless p > P, in which 
case we proceed to step 9. 

9. Let r be the least residue in absolute value of E (mod A). 

10. Return the value of r, and STOP. 

It is relatively easy to see that the algorithm requires no more than the 
indicated space and time requirements. 

The idea behind the algorithm is that we use the Chinese Remainder 
Theorem to compute the x''^^ coefficient of our polynomial, and the compu- 
tation in step 7 is just an "on the fly" CRT calculation. This calculation is 
based on the following fact: If qi, ...,qh are coprime, and if ai, ...,ah are any 
integers, then if we set 

i=l 

and 

A' 

T! = — , wYieie hi = ai{A' / qi)~^ (mod gj), 

then 

E' = Oj (mod qi), for every i — 1,2, h. 

One might guess that the coefficient of our polynomial can be computed 
using less resources by using a "Fourier Series" method; that is, 

Coei. x"-' in fix) = i |j e-''^^^M/^/(e2'^'^/^). 

£=0 

It is not obvious (to me) how to do this without using the special form of 
the polynomial f{x): First of all, we would need to maintain ^ T digits 
of precision for each term in the sum, since any particular coefficient of the 
polynomial f{x) can have size 2'^^^°^-^, for some c > 0. Thus, 3> bit 
operations would be needed to compute each term /(e^'^'^/^). In total, » 
bit operations would be needed to evaluate all the terms in the sum. If one 
tries to use FFT's to evaluate all the terms in the sum at the same time, 
this reduces the running time to O(T^log^T) bit operations; however, the 
memory requirements then increase to 3> bits of storage, which is the 
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amount needed to store the all numbers /(e^'"*^/"^), 0<i<Tto^T bits 
of precision. Even if we try a discrete version of this method, where the 
polynomials are computed, say, modulo 2^ for k ^ T, and the roots of unity 
are roots of unity modulo 2'^, we would run into the same difficulties. 



4 Proof of Proposition [T] 

The Aj's can be computed by solving the equation 



/ \ 



As 



\AtJ 

If we call the matrix on the left-hand-side M, then 

A, = na^-M-l 



n 


1 


1 


1 


\ 


1 


2 


3 


T 




1 




32 . 


■■ T2 




Vi 


2T-1 


3T-1 . 




V 





no!' 




(9) 



V / 



(10) 



where I is the entry in the jth row, fcth column of M ^ . 

We will calculate M~l via polynomial interpolation: We have that for 
any set of ordered pairs 

(1,61), (2,62), (T,M, 

where 61, &t ^ C, there exists a unique degree T— 1 polynomial f\x) G C[x] 
such that 

f{i) = bi, for alH = 1,2, ...,T; 

moreover, if we write 



f{x) = CtX + Ct-iX H h C2X + Ci 
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then these coefficients c, can be calculated in two different ways: The ffist 
way is through basic linear algebra, since 



(11) 



We notice that the matrix on the left-hand-side is M', the transpose of our 
matrix M. 

The second way of calculating the Cj's is by Lagrange interpolation, which 
gives 

T T 



(l 


1 


1 ■ 


1 


\ 


/Ci\ 






1 


2 


22 . 


■■ 2^' 


-1 


C2 






1 


3 


32 . 


■■ 3^- 


-1 


C3 






Vi 


T 


T2 . 


.. T^- 


v 


Vt) 







Now, if we suppose that 



i=l h^l 



X — h 
i — h 



0, ifz^j, 
_ 1, if ^ = j, 

then for this choice of 6j's, one sees from (II Ij) that 

Ck = (M%-J = M-^ 
On the other hand, from (|12j) we see that 

T 



Ck = Coef. of x'' Mn JJ^ 



X — h 
j — h 



Thus, 



Coef. of a;'' Mn JJ^ 



X — h 



and we conclude from this and (fTUI) that ((Tj) holds. 
Finally, to prove (jHl), we note that 



h = l 



(12) 
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The coefficient of x'^ ^ in the above polynomial is clearly less than 



which proves (jH))- 
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